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We present a new Skyrme-Hartree-Fock-Bogoliubov nuclear-mass model in which the contact 
pairing force is constructed from microscopic pairing gaps of symmetric nuclear matter and neu- 
tron matter calculated from realistic two- and three-body forces, with medium-polarization effects 
included. With the pairing being treated more realistically than in any of our earlier models, the 
rms deviation with respect to essentially all the available mass data falls to 0.581 MeV, the best 
value ever found within the mean-held framework. Since our Skyrme force is also constrained by 
the properties of pure neutron matter this new model is particularly well-suited for application to 
astrophysical problems involving a neutron-rich environment, such as the elucidation of the r-process 
of nucleosynthesis, and the description of supernova cores and neutron-star crusts. 
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With a view to their astrophysical application in 
neutron-rich environments, we have developed a se- 
ries of nuclear-mass models based on the Hartree-Fock- 
Bogoliubov (HFB) method with Skyrme and contact- 
pairing forces, together with phcnomcnological Wigner 
terms and correction terms for the spurious collective en- 
ergy. All the model parameters are fitted to essentially 
all the experimental mass data. 

Model HFB-9 [ij and all later models constrained the 
underlying Skyrme force to fit the energy-density curve 
of neutron rnatter, as calculated by Friedman and Pand- 
haripande [3] for realistic two- and three-nucleon forces. 
In the latest of our published models, HFB-16 we 
imposed a comparable constraint on the contact pairing 
force. Instead of postulating a simple functional form for 
its density dependence, as is usually done, we constructed 
the pairing force by solving the HFB equations in uniform 
matter and requiring that the resulting gap reproduce ex- 
actly, as a function of density, the microscopic ^5*0 pairing 
gap calculated with realistic forces. In that preliminary 
study we assumed that the pairing strength for neutrons 
(protons) depended only on the neutron (proton) density, 
as suggested by Duguet [3] , and chose for this microscopic 
reference gap the one calculated for pure neutron matter 
without medium effects Q. We obtained thereby what 
was at the time our best-ever fit to the mass data, the 
rms deviation for our usual data set of 2149 measured 
masses of nuclei with N and 2' > 8 [6] being 0.632 MeV. 
On the other hand, the mass fits were much worse if we 
chose reference pairing gaps calculated with medium ef- 
fects taken into account. 

Here we show that it is possible to obtain excellent 
mass fits even when the pairing force is constrained to 
microscopically calculated gaps in which medium effects 
have been included. The essential step is to impose the 
additional constraint of asymmetric nuclear matter pair- 



ing, thereby allowing the neutron and proton pairing 
strengths each to depend on both the neutron and proton 
densities. 

The HFB- 17 mass model. With this generalization of 
our earlier pairing model, Eq. (3.3) of Ref. 0] for the 
pairing strength is replaced by 

/ \ 

V^'iPn, Pp] = -StT^ — — X 

\2M*{pn,Pp)J 

where Aq{pn, pp) is the corresponding pairing gap of 
asymmetric nuclear matter calculated microscopically, 
M*{pn, Pp) is the effective nucleon mass and e\ is the 
pairing cutoff. The chemical potential pq is approxi- 
mated by pq = ri^klJ{2M*), where fcpq = [iir'^Pqf/^ 
is the Fermi wave number. 

For the reference microscopic gap we use the recent 
Brueckner calculations of Cao et al. which are based 
on realistic two- and three-nucleon forces. All these cal- 
culations include the effect on the interaction of medium 
polarization, and are performed both with and without 
self-energy corrections. Ideally, we should have used the 
former gaps, which imply an effective mass M* different 
from M . However, for complete consistency in Eq. ([T|) 
our Skyrme force would then have been required to re- 
produce the same M*. Satisfying this further constraint 
in addition to all the others that we have already imposed 
seems to be impossible within the framework of the con- 
ventional Skyrme forces used here. Thus if one wishes to 
retain the self-energy corrections the best that one can 
do is to relax this constraint on the effective mass, but 
we found that this option was incompatible with good 
mass fits. We thus take the gaps calculated without self- 
energy effects, and then for consistency set M* = M in 
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FIG. 1: (Color online) ^So pairing gap A in infinite neutron 
(square) and symmetric (circle) nuclear matter as a function 
of the Fermi wave number, taken from Ref. p]. 



Eq. II]). This was the choice that led to the exceUent 
mass fits reported below. 

Ref. [3] calculates pairing gaps only for symmetric nu- 
clear matter, As nip = Pn + Pp), and pure neutron mat- 
ter, ANAiipn)- Since we need the pairing gaps for arbi- 
trary asymmetry we adopted the interpolation ansatz 



\iPn,Pp) = Asm(p)(1 - |?7l) ± ANM{Pq)v 



Pq 



(2) 



where 7y = (p„ ^ Pp)l P and the upper (lower) sign is to 
be taken for q = n{p)] we have also assumed charge sym- 
metry, i.e., An{pn,Pp) = Ap{pp,pn). This expression en- 
sures that for symmetric nuclear matter, Aq{p/2, p/2) = 
^Sm{p) and for neutron matter A„(/9, 0) — Anm{p) and 
Ap(p,0) = 0. 

Because of Coulomb effects, and a possible charge- 
symmetry breaking of nuclear forces, we must allow for 
the proton pairing strength to be different from the neu- 
tron pairing strength. Likewise, we follow our usual prac- 
tice of allowing the pairing to be slightly stronger for nu- 
cleons of which there are an odd number. This procedure 
can be understood microscopically Q as compensating 
for the neglect of the time-odd fields implicit in our use of 
the equal- filling approximation (EFA) [9[ . (Note that the 
odd nucleon will nevertheless contribute to the time-even 
fields.) We take account of these extra degrees of free- 
dom by multiplying the value of v'^'^[pq], as determined 
through Eq. ([1]), by renormalizing factors where 
fp,fp and /~ are free, density-independent parameters 
to be included in the mass fit. We set /+ = 1, tac- 
itly supposing that all effects related to charge-symmetry 
breaking act only on protons. 

Results. The foregoing model, labeled HFB-17, was 
fitted to the above data set of 2149 measured nuclear 
masses @; in making this fit we followed our recently 
adopted strategy of dropping the Coulomb-exchange 
term. The resulting parameter set, labeled BSkl7, is 
given in Table U (definitions of these parameters can be 
found in Q). 
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FIG. 2: Differences between measured |6|] and HFB-17 masses, 
as function of A'^. 



TABLE L Parameter set BSkl7: left panel shows the Skyrme 
parameters, middle panel the pairing parameters and right 
panel the parameters for the Wigner and collective corrections 
(units for energy and length are MeV and fm respectively). 
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The deviations (data-theory) between all the 2149 
measured masses of our data set and the new predictions 
are shown graphically in Fig. [21 no deviation exceeded 2.8 
MeV. The rms and mean values of these deviations are 
shown in the first two lines of Table HIl with an rms devi- 
ation of 0.581 MeV this is the most accurate mass model 
ever achieved within the mean-field framework. The next 
six lines of this table refer to various subsets of our data 
set. We stress that all the 2149 data points to which we 
make our fit are taken from the 2003 Atomic Mass Eval- 
uation (AME) However, a considerable amount of 
mass data has accumulated in the meantime, but since 
these new measurements have not been subjected to the 
same scrutiny that went into the 2003 AME we have ex- 
cluded them from our fit. Nevertheless, it is of interest to 
compare these new data with our model, and we do this 
in lines 9 to 12 of Table [TTl for two sets of measurements, 
Refs. Uj and 1^. It is remarkable that our model agrees 
better with these new data than with the fitting data. 

We have pointed out in Ref. lj| that only slight 



changes in the quality of the mass fit can be accompanied 
by considerable changes in the overall pairing strength, 
as measure d by the spectral pairing gaps (defined in Eq. 



(5) of Ref. \XM)- Now throughout our entire mass-model 



project we have been concerned not only with nuclear 
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TABLE II: Rms (a) and mean (e) deviations between data [g] 
and HFB-17 predictions. The first pair of lines refers to all 
the 2149 measured masses M, the second pair to the masses 
Mnr of the subset of 185 neutron-rich nuclei with S„ < 5.0 
MeV, the third pair to the 1988 measured neutron separation 
energies Sn and the fourth pair to 1868 measured beta-decay 
energies Q/3. The fifth and six pairs correspond to the devia- 
tion with respect to the recently measured masses of Ref. [ill ] 
and [l2|, respectively. The seventh pair shows the compari- 
son with the 782 measured charge radii [l3|, and the last line 
shows the calculated neutron-skin thickness of ■^°*Pb. Note 
that units for energy and length are MeV and fm respectively. 





HFB-16 HFB-17 


cr(2149 M) [6] 


0.632 
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e(2149 M) [6] 
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masses but also with nuclear level densities (among other 
quantities of astrophysical interest), and these are very 
sensitive to the spectral pairing gaps. The values of these 
quantities found here are comparable to those of model 
HFB-13 IJI, which lead to satisfactory level densities. 

We have determined the ground-state spins and pari- 
ties of odd nuclei using the HFB-17 single-particle level 
scheme obtained in the EFA, as described above. For 
odd-A nuclei, the spin and parity are assumed to be those 
of the single nucleon of the last filled orbit, while for odd- 
odd nuclei we apply the Nordheim rule [15]. For the 90 
spherical odd-A nuclei with quadrupole deformation pa- 
rameter /32 < 0.05, 91% of the experimental spins [16| 
are correctly predicted, while for the 717 deformed ones 
with > 0.16, only 41% are correctly determined. For 
all the 1582 odd-A and odd-odd nuclei we obtain a global 
success rate of 47% on the spins and 72% on the parities 
(we assume a spherical configuration for /32 < 0.16). The 
most important ground state properties predicted by the 
HFB-17 model, including spins and parities, have been 
tabulated for all the 8508 nuclei with 8 < Z < 110 be- 
tween the proton and neutron drip lines. 

Table HH] shows the macroscopic parameters (infinite 
and semi-infinite nuclear matter) calculated for the force 
BSkl7 (for the definition of these parameters see, for ex- 
ample, Ref. Q). This table also shows the values of these 
parameters for force BSkl6, the force underlying mass 



TABLE III: Macroscopic parameters for forces BSkl6 and 
BSkl7. The first twelve lines refer to infinite nuclear matter, 
the last two to semi-infinite nuclear matter. Note that units 
for energy and length are MeV and fm respectively. 
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FIG. 3: (Color online) Energy per neutron (MeV) as a func- 
tion of density (fm~'^) of neutron matter for BSkl7 and for 
the calculations of Ref. |^ (FP) and A18-|-^i;-|-UIX* of Akmal 
et al. Mi. 



model HFB-16, and it will be seen that in this respect 
there is very little difference between the two forces. This 
is hardly surprising given that the macroscopic parame- 
ters depend entirely on the Skyrme force, and it is in the 
pairing channel that we have introduced the principal 
modifications. (Note that the values of the isoscalar ef- 
fective mass M* and symmetry energy J were imposed.) 
It will be seen that in both models the isovector effective 
mass M* is found to be smaller than M* at the saturation 
density po, implying thereby that the neutron effective 
mass M* is larger than the proton effective mass M* in 
neutron-rich matter. Such an isovector splitting of the ef- 
fective mass is consistent with measurements of isovector 
giant resonances 17|, and has been confirmed in several 



many-body calculations with realistic forces [18j, [19| . 

As seen in Fig. [31 the energy-density curve of neutron 
matter for force BSkl7 is identical to the realistic curve of 
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Ref. [3| up to the supernuclear density of 0.4 fm ^, as is 
the case with ah our other forces that have been fitted to 



J = 30 MeV. It is to be noted that unhke Ref. [17| we have 
not had to resort to a second term in the Skyrme force 
in order to simultaneously fit neutron matter and obtain 
the correct sign for the isovector splitting of the effective 
mass. Fig. [3] also shows the energy- density curve given 
by the realistic calculation AlS+^u+UIX* of Akmal et 
al. [iSl- There is more physics in this calculation than in 
the one of FP, but we are unable to fit this curve without 
degrading the quality of the mass fit. However, there 
have been some very recent indications that this curve 
might be too steep l2l|. 

Fig.[4]shows the potential energy per particle in each of 
the four two-body spin-isospin {S, T) channels as a func- 
tion of density for symmetric nuclear matter; we give re- 
sults for both BSkl7 and Brueckner-Hartree-Fock (BHF) 
calculations with realistic two- and three-nucleon forces 
[i^ l- A fair agreement between BSkl7 and the realis- 
tic calculations in all states can be seen; note particu- 
larly that the deviation in the (1,1) channel is much less 
marked than in Ref. [s*!, mainly because there we com- 
pared with older BHF calculations. 

Conclusions. We have described a new Skyrme-HFB 
nuclear-mass model, HFB-17, in which the contact pair- 
ing force is constructed from microscopic pairing gaps of 
symmetric nuclear matter and neutron matter calculated 
from realistic two- and three-body forces, with medium- 
polarization effects included. In this way the rms de- 
viation with respect to essentially all the available mass 
data has been reduced, for the first time with a mean-field 
model, below 0.6 MeV. Given also the constraint imposed 
on the Skyrme force by microscopic calculations of neu- 
tron matter, this new model is particularly well adapted 
to astrophysical applications involving a neutron-rich en- 
vironment, such as the elucidation of the r- process of nu- 
cleosynthesis, and the description of supernova cores and 
neutron-star crusts. 
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